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In all iron pnictides, the positions of the ligand alternatively above and below the Fe plane create 
2 inequivalent Fe sites. This results in 10 Fe 3d bands in the electronic structure. However, they do 
not all have the same status for an ARPES experiment. There are interference effects between the 2 
Fe that modulate strongly the intensity of the bands and that can even switch their parity. We give a 
simple description of these effects, notably showing that ARPES polarization selection rules in these 
systems cannot be applied by reference to a single Fe ion. We show that ARPES data for the electron 
pockets in Ba(Feo.92Coo.os)2As2 are in excellent agreement with this model. We observe both the 
total suppression of some bands and the parity switching of some other bands. Once these effects are 
properly taken into account, the structure of the electron pockets, as measured by ARPES, becomes 
very clear and simple. By combining ARPES measurements in different experimental configurations, 
we clearly isolate each band forming one of the electron pockets. We identify a deep electron band 
along one ellipse axis with the d xy orbital and a shallow electron band along the perpendicular axis 
with the d xz /dy Z orbitals, in good agreement with band structure calculations. We show that the 
electron pockets are warped as a function of k z as expected theoretically, but that they are much 
smaller than predicted by the calculation. 

PACS numbers: 79.60.-i, 71.18.-y, 71.30.-h 



I. INTRODUCTION 

There is a consensus that the multiband nature of 
iron pnictides is essential for defining their electronic 
properties. Many models for the superconducting and 
antiferromagnetic orders heavily rely on the interaction 
between different electron and hole Fermi Surface (FS) 
sheets^! However, despite many ARPES investigations 
of iron pnictides, the structure of the electron pockets 
are still poorly understood. They appear more diffi- 
cult to clearly detect in ARPES than the hole pockets. 
Usually, only part of the expected ellipse is observed 3 -"^, 
leaving unclear what is really the size of the pocket or 
whether the electronic properties are isotropic around 
the ellipse. Recently, an ARPES study resolved them in 
more details 6 -', but proposed a structure so different from 
band structure calculations, that it clearly calls for more 
investigation. On the contrary, other experimental tech- 
niques, such as trans port!^ , RamaiP^ or quantum os- 
cillations experiment d 11 * 12 !, seem to detect predominantly 
electrons, as if they had longer lifetimes. It would be de- 
sirable to measure the Fermi velocities and lifetimes of 
each hole and electron bands separately by ARPES, to 
understand the origin of this difference. This has not 
been possible so far, due to the low resolution obtained 
on electron pockets. 

In this paper, we present a complete investigation of 



the electron pockets in Ba(Feo.92Coo.os)2As2. This com- 
pound corresponds to optimal Co doping for supercon- 
ductivity (T C =23K), but we have obtained similar spec- 
tra for BaFe2 Assr^ and other Co dopings. A complexity 
of ARPES in iron pnictides is that there are 2 Fe per 
unit cell, which doubles the number of bands. In section 
II, we give a simple description of the folding from the 
lFe to 2Fe Brillouin Zones (BZ), which yields a quali- 
tative understanding of the spectral weight distribution 
on the main and folded bands. To our knowledge, fold- 
ing effects were only discu ssed b efore in ARPES in the 
context of a weak potential 16 1 17 1 . This case is different, 
because the As potential is strong, but essentially only 
its symmetry changes between the 2 inequivalent Fe sites. 
We show that this does modulate strongly the spectral 
weight, but in a way that depends on the orbital sym- 
metry. Understanding this allows to identify the bands 
that should be observed in ARPES. This also explains 
that the parities of the bands depends on the relative 
phase between the 2 Fe of the unit cell. This implies that 
the polarization selection rules for ARPES are consider- 
ably different from those generally used in this commu- 
nity (see section II-C) and some bands appear as if they 
have "switched" parity. In section III, we present our 
ARPES data, where one full electron ellipse is clearly re- 
solved. We first compare our measurements with band 
structure calculations, which allows to identify the band 
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Figure 1. (a) Structure of one FeAs slab. The thick black line represents the 2D unit cell with inequivalent Fe at the center 
(A) and the corners (B), because of the position of the neighboring As. The dotted black line defines the basic Fe square, (b) 
Tetrahedral As environment of one Fe. (c) 2D Fermi Surface expected in the 2Fe BZ built on the unit cell. The color code for 
different orbital characters. The pockets at F and M are hole pockets, those at X are electron pockets. The dotted square is the 
lFe BZ built on the basic Fe square, (d) FS in the lFe BZ. The additional pockets in the 2Fe BZ can be obtained by folding 
with respect to the 2Fe BZ boundaries or, equivalently, by translation by qF- We follow here ref. 13 for the definition of main 
and folded pockets, other authors reverse the T and M points^, (e) Sketch of the expected spectral distribution, following the 
analysis of section II. 



characters. Using the correct polarization selection rules 
solves the problems encountered in ref. 6 and reconciles 
ARPES with band calculations. We then compare in sec- 
tion III-B the measurements with simulations using the 
unfolded band structure of refffl This gives the spectral 
weights of the main and folded bands, which are found 
in excellent agreement with our measurement. We finally 
study in III-C the k z dependence of the different electron 
bands. The size of the pockets and their k z dependence 
are smaller than predi cted in theory, confirming a ten- 
dency previously notecP^^Sl. A shift of the calculated 
electron bands as large as 150meV is needed to fit the 
experiment. 



Single crystals were grown using a FeAs self-flux 
method and were studied in details by transport 
measurements 7 . ARPES experiments were carried out at 
the CASSIOPEE beamline at the SOLEIL synchrotron, 
with a Scienta R4000 analyser, an angular resolution 
of 0.2° and an energy resolution better than 15 meV. 
Band structure calculations were perfomed within the lo- 
cal density approximation, using the Wien2K package^, 
with the experimental structure of BaFe2As2 and a dop- 
ing of 10% treated in the virtual crystal approximation. 



II. BANDS IN THE 1 FE AND 2 FE BZ 

In all iron pnictides, there are 2 inequivalent Fe per 
unit cell, as sketched in red and blue in Fig. 1(a), and 
called sites A and B in the following. The inequivalency 
is due to the position of the surrounding As, which always 



form a tetrahedra, but with the edge oriented along 110 
either above the Fe plane, as in Fig. 1(b), or below. The 
As potential created at the 2 inequivalent Fe sites is equal, 
as first approximation, but with an opposite symmetry 
with respect to z. 

In the 2Fe BZ corresponding to this structural unit cell 
[thick line in Fig. 1(c)], there are 10 Fe 3d bands to con- 
sider, making the band structure quite complicated to 
decipher. To simplify this problem, many theoreticians 
use a lFe BZ, as in Fig. 1(d), which is based on the small 
Fe square [dotted line in Fig. 1(a)], and where the num- 
ber of bands is divided by two. It is well known that the 
other bands can be obtained to a good approxim ation 
by folding with respect to the 2Fe BZ boundarie J 1 | 13 | 14 [ 
Curiously, very little attention has been paid to this un- 
derlying structure by ARPES authors. It is usually con- 
sidered as a simple convention to work in the lFe or 2Fe 
BZ. On the contrary, the distribution of spectral weight 
on the main (lFe) and folded bands is an intrinsic fea- 
ture of the electronic structure. We will show that this 
modulates strongly the ARPES intensities of the bands 
and that it is very useful to understand how this works 
to correctly interpret ARPES spectra. 

In the following, we choose the x and y axis in the Fe- 
Fe direction, as in Fig. 1(a), the origin at one Fe site, 
and we call a the Fe-Fe distance. The orbitals are also 
defined with respect to these axes, which are those of the 
lFe BZ. In Fig. 1, we consider only 2 dimensional (2D) 
BZ. In reality, there is a component of the folding along 
k z for BaFe2As2, because it is body centered^. This 
will be considered later in section III. We label the high 
symmetry points X and M according to the 2Fe BZ, to 
keep the same notations between sections II and III. 



A. Folding bands from the lFe to 2Fe BZ 

We first consider a square lattice of Fe. We define 
Bloch functions ^\ F ^(r) built from atomic x( r ) orbitals. 

^l F x e (r) = J2e ik - R x(r-R) (1) 

R 

x(r) is one of the 3d Fe orbitals and the sum runs over 
all Fe sites (for clarity, we will omit the index \ in the 
following). At T, all orbitals are in phase and a spatial 
representation of ^\ Fe {r) for d xy is given in Fig. 2(b). 
At X, i.e. k=(7r/a,0), there is an additional dephasing 
(j) = k ■ R due to the Bloch term. For the Fe sites along 
y, (j) = and for the other sites at (±a,0), <j) = tt. This is 
sketched in Fig. 2(c), again for d xy . 

The BZ of this square lattice is the lFe BZ. How- 
ever, one could decide to work in the 2Fe BZ and obtain 
an equivalent description of these bands by folding ifj 1Fe 
with respect to the new zone boundaries, or, equivalently, 
translating them by q — (ir/a,ir/d) [see Fig. 1(d)]. This 
creates a new band ^°^ erf ( r ). 

^oi ded{r) = ^iFe {r) = J2 e ^ R e ik - R X {r - R) 

R 

For the Fe on sites A, R = (0,0), yielding e iq,R = 1. 
For sites B, R = (±a, 0) or (0,±a), yielding e iq ' R = -1. 
Hence, 

^ Folded (r) = J2e ik ' RA x(r - R A ) - J^^'^xir " r b) 

Ra Rb 

In other words, the Fe on sites B are dephased by tt 
compared to sites A. This is represented in Fig. 2(d) 
for d xy at T. An equivalent description of the 2 bands in 
the 2Fe BZ is then to define a linear combination of the 2 
Fe with sites A and B in-phase or out-of-phase. We will 
call ip the in-phase bands and -0* the out-of-phase bands, 
like d xy in Fig. 2(b) and d xy * in Fig. 2(d). 

The choice of the lFe or 2Fe BZ is of course dictated 
by physical considerations. Usually, the 2Fe BZ is re- 
quired by the presence of a small potential V making the 
2 Fe in the unit cell inequivalent. The eigenstates ip 1 
and ip 11 will be a linear combination of ip 1Fe and ^ Folded 
(or equivalently of ip and The coupling will for ex- 
ample open gaps at the bands crossings, proportional to 
the folding potential V. In iron pnictides, the real bands 
may also hybridize with other orbitals having the same 
symmetry and acquire some finite weight from different 
orbital character. Nevertheless, if these perturbations 
remain sufficiently small, the eigenstates behave essen- 
tially as the bands we have defined here, and they keep 
the same symmetry, so that these bands are very useful 
guides to simply visualize the situation. 

B. Spectral weight distribution 

Since the final state of the photoelectron of the ARPES 
process corresponds to a free-electron like state far away 
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Figure 2. (a) Sketch of the three t2 9 orbitals, d xy , d xz and 
d yz - (b-e) Relative orientations of orbitals in real space for 
the following orbitals and points of the reciprocal space, (b) 
d xy at T, (c) d xy * at X(7r/a,0), (d) d xy * at T, (e) d xz * at X. 
The blue line in (b) indicates the xz mirror plane. 



from the sample, the ARPES spectral weight is natu- 
rally given by the lFe picture with similar translational 
symmetry. The spectral weights w 1 ^ 11 of the 2 eigenstates 
^ / ' // can be obtained by the projection on the basis of 
the lFe BZ. With the previous notations, 

<' U = \<ri Fe \ri n >\ 2 (2) 

For a weak folding potential, the mixing between ip 1Fe 
and ifj Folded is rather small, except near band crossings, 
and the spectral weight remains concentrated along the 
bands of the lFe BZ. Many such examples can be found 
m ARPES literature™!. In iron pnictides, the potential 
due to the As is certainly not weak, as the As is essential 
to define the overlap between the Fe^. On the other 
hand, this potential is very similar at each Fe site, in 
fact, it is essentially reversed with respect to z. This 
leads to original folding effects, where half the Fe bands 
do have weaker spectral weight, but not those defined as 
the lFe bands in Fig. 1(d). 

To perform the calculation of Fig. 1(d), all Fe must 
be equivalent. As the As potential cannot be neglected, 
a way to achieve this is to introduce a virtual operation 
symmetry, which changes the sign of z when going from 
Fe A to B. This is the glide- mirror symmetry described by 
Andersen and BoerP^ under which all Fe become equiv- 
alent. As d xy is symmetric with respect to z [see Fig. 
2(a)], it does not change between the two sites under 
this glide-mirror symmetry. On the contrary, d xz and 
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Figure 3. (a) Sketch of the bands forming the electron pock- 
ets. Thick and thin lines indicate main and folded bands in 
the sense of the lFe BZ with glide-mirror symmetry. Solid 
and dotted lines represent the in-phase and out-of-phase or- 
bitals. Colors indicate the main orbital character, (b) Sketch 
of the dispersion of the different bands along the path TXM. 
The color sketches the main orbital character, although there 
may be hybridization with other orbitals. (c) Symmetries of 
the different bands forming the electron pocket at X. 



d yz change sign. It follows that the "main bands" of this 
lFe BZ are the one we defined as in-phase for d xy , but 
out-of-phase for d xz and d yz . Eq. (2) predicts that the 
spectral weight will be close to 1 for all in-phase bands 
and to for all out-of-phase bands. This will give rise to 
incomplete ellipses and inequivalency between T and M, 
as sketched in Fig. 1(e). This is also a lFe representation 
of the electronic structure, but where the real space co- 
ordinate has not been redefined by utilizing glide-mirror 
symmetry. Instead, a flat reciprocal space of lFe BZ has 
been constructed by treating the translational symmetry 
breaking as the perturbation. This definition retains the 
symmetry of the physical space where ARPES works. 

The recently proposed band structure "unfolding" 
methocP^l provides a simple mean to calculate such an 
unfolded band structure from first-principles calculation. 
By representing the one-particle spectral function in the 
lFe Bloch basis, the above phase interference effect (as- 
sociated with the Fe positions) that was hidden in the 
2Fe Bloch functions now enters explicitly the theoretical 
spectral weight. This was recently worked out explicitly 
in ref. EH. Essentially, it finds that the spectral weight 
is concentrated along the in-phase bands, in agreement 
with our qualitative reasoning. In addition, it gives the 
residual intensity for the out-of-phase bands, which is 
found to strongly depend on the hybridization with other 
orbitals. Indeed, we have shown that bands symmetric 
with respect to z will give weight along the in-phase band 
and vice- versa. If two bands of different symmetry with 
respect to z are hybridized, their weight will appear on 
different types of bands (in-phase or out-of-phase). For 
example, the d xy band is often significantly hybridized 
with d xz /d yzi and the out-of-phase bands then have sig- 



nificant weight with mainly d xz /d yz character. This was 
called "parity switching" in refffl. 



C. Band symmetry 

One important consequence of defining the wave func- 
tions the way we have done is to be able to determine 
their symmetry. Along k x , the important symmetry op- 
eration is the parity with respect to the xz mirror plane. 
In Fig. 2(b), one can notice that the xz plane containing 
the Fe is not a mirror plane (the As are not symmet- 
ric with respect to that plane). On the other hand, the 
one containing the As (blue line) is a true mirror plane. 
Therefore, the symmetry of the wave functions does not 
only involve the orbital character, but also the relative 
phase between the 2 inequivalent Fe. This situation is 
really a direct consequence of the existence of 2 inequiva- 
lent Fe that is induced by the As positions, as the mirror 
plane would contain the Fe, if there was no As. In-phase 
and out-of-phase bands then have opposite symmetries 
by construction, as sketched in Fig. 2(b) and 2(d) for 
d xy and d xy *. This "parity switching" of d xy * also defines 
the symmetry of the orbitals with which it can hybridize. 
As we have seen that the d xy weight will be suppressed 
by unfolding effects, these hybridized orbitals will be re- 
sponsible for the residual weight. 

The symmetry also changes with the position in recip- 
rocal space, as the Bloch phase term changes. This is 
illustrated in Fig. 2(c) and 2(e) for d xy and d xz * at X, 
d xy is for example odd/yz at Y but even/yz at X. Here 
again, if the Fe were in the mirror plane, the change 
in phase in k-space would not change the band symme- 
try, so that it is a consequence of the 2 Fe situation. 
In Fig. [3ja), we sketch the electron pockets at X, coded 
with all features that will be important to understand 
the spectral weight distribution and symmetries : the 
main (thick lines) and folded (thin lines) characters, in 
the sense of the lFe BZ with glide-mirror symmetry [Fig. 
1(c)]; the in-phase (solid lines) and out-of-phase (dotted 
lines) characters, as defined in II-A. The dispersion and 
symmetries with respect to the main mirror planes xz 
and yz are also indicated in Fig. [3jb) and (c). To deter- 
mine if one band is in-phase or out-of-phase, we use Fig. 
1(d) where the represented d xy (resp. d xz /d yz ) bands 
are in-phase (resp. out-of-phase), as explained before. 
The dispersions are illustrated by the band structure of 
BaFe2As2, but the symmetry features remain the same 
for all structures, only the k z dependences are different. 
One can see that the 2 bands forming the folded electron 
pockets, in the sense of the lFe BZ with glide- mirror sym- 
metry, d yz and d xy *, have opposite symmetries along TX 
and the same symmetries along XM. Consequently, they 
cross without hybridizing along TX, but repel each other 
along XM. This can only be understood after considering 
the symmetries derived here. As ARPES selection rules 
depend on the parity with respect to the mirror planes, 
it is also crucial to understand this to work this out to 
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apply them correctly. 



III. ARPES ON ELECTRON POCKETS 

In Fig. [4j we now use the true structure of BaFe2As2, 
whose 3D BZ is sketched in Fig.[4^c). The only difference 
is that the M point in Fig. 1 is now equivalent to Z and 
q connects different k z . In Fig. Eta), the FS is sketched 
for k z =l, where the FS of Fig. ]4[b) was measured (see 
below). Because of this, the two ellipses are now ori- 
ented in the sa me d irection, instead of being perpendic- 
ular as in Fig. J3[ 4 . Only the k z dispersion make their 
shape slightly different (they are slightly more rounded 
at /c 2 =CM^, more information on k z dependence will be 
given in section III-C). 



A. Assignment of the different bands 

Fig. [4^b) presents the FS at 25K, at a photon 
energy fkj=34eV. This corresponds at normal emis- 
sion to k z =l modulo tt/c\ where c'=6.5A is the dis- 
tance between 2 FeAs planes. Indeed, the photon 
energy Hlj selects one particular k zi through k z = 

yj2m/h 2 * (hoj — W + Vb) — Ay/, where m is the electron 

mass, kj i the momentum in xy plane, W the work func- 
tion and Vb an inner potential, fixed here at 14eV, as 
already commonly used for pnictide d 4 l 20 i 

Two circular hole pockets are detected at Z, as usually 
observed in the literature. In most ARPES studies pub- 
lished so faiP^, the electron pockets look like in panel 
a, with the sides of the expected ellipse oriented towards 
Z missing. In panel /?, however, we resolve clearly the 
parts towards Z, while the middle parts are missing. In 
Fig. [4je), we combine the FS crossings detected in these 
two panels (green points for a and blue points for j3. The 
latter are rotated by 90°, as k x and k y are equivalent by 
symmetry). They define a nearly perfect ellipse, match- 
ing the one sketched by thin line in Fig. [4^a). 

The difference between a and /3 is due to the direction 
of light polarization, which defines the symmetry of the 
orbitals that can be observ ed. A s already widely used 
in ARPES of iron pnictide d 3 * 6 * 24 *, to select even or odd 
orbitals with respect to a mirror plane of the structure, 
the polarization must be even or odd with respect to 
that plane 25 . In our experiment, the polarization is fixed 
along k y and is therefore odd/xz and even/yz [see Fig. 
|4^d)]. To view odd orbitals, we detect electrons in the 
xz plane and, for even orbitals, we detect electrons in yz. 
This type of geometry was for example used in refpa 
although it is more common to detect electrons in one 
plane and switch the polarization from in-plane (so-called 
p configuration, selecting even orbitals) to out-of-plane (s 
configuration for odd orbitals). For our case, the inten- 
sity of the even band is much lower in p-polarization. 
If Fe would belong to the mirror plane, these selection 



rules would only depend on the orbital character. Such a 
situation was implicitly assumed in all previous studies. 
As discussed previously, this is in fact not the case here 
and one should use the correct symmetries determined 
in Fig. [3jc) , which change both as a function of location 
in the reciprocal space and of the in-phase/out-of-phase 
character. It turns out that all hole pockets at T (or Z 
along 001) have in-phase characters, so that the selection 
rules just follow the orbital symmetry for hole pockets. 
However, they are completely different for electron pock- 
ets. According to the band symmetries defined in Fig. 
Kmc), the allowed orbitals in a are d xz *, d xy and d yz . For 
p, these are d xz *, d yz and d xy *. If we rotate this panel 
again by 90° to compare it directly to a, we just have to 
switch the roles of x and y and we expect to see orbitals 
even/xz, which are d yz *, d xz and d xy *. 

In Fig. [5j we show the dispersions of the bands along 
the cuts indicated in Fig. |4je), the major ellipse axis in 
(a), the minor ellipse axis in (b) and an axis parallel to the 
major axis at k y =0.0b A -1 in (c). For even bands, we use 
again panel f3 rotated by 90°. We observe that the even 
electron band along the major axis [Fig.l5jal)] is very dif- 
ferent from that along the minor axis [Fig. f5^b2)] . It is 
much deeper (at least lOOmeV) and with steeper disper- 
sion, while the odd band is quite shallow (50meV). The 
calculation shows that this large anisotropy is intrinsic to 
the orbital character. The shallow band is mainly formed 
by d xz or d yz orbitals and the deep band by the d xy or- 
bital. This difference in dispersion makes it very easy 
to pick up the band character; it is d yz for the shallow 
band and d xy * for the deep band. d yz is odd/xz, so that 
we expect to see it in a, as we do. d xy * is even/xz and 
appears in /3, as we expect. Note that this symmetry is 
opposite to that of the d xy orbital, so that this complete 
analysis is necessary to understand its parity. Because of 
this problem, it was wrongly assigned to d X 2- y 2 in refP. 
In ref. 15 , we correctly assigned it to d xy , because of its 
dispersion, but we did not understand its parity. 

Qualitatively, we predicted in section II the spectral 
weight to be concentrated along in-phase bands (solid 
lines). This explains very well the strong weight of d yz in 
(a2-c2) and also the strong suppression of d xz *. In fact, 
this suppression is direct evidence for the crucial impor- 
tance of these unfolding effects. We would also expect to 
observe clearly d xy , for example in (b2), which is how- 
ever unclear. We attribute the absence of d xy to its much 
smaller cross section for ARPES. Indeed, the smallness of 
polar angle of the emitted photoelectron suppresses the 
dipole matrix element magnitude of the purely in-plane 
orbitals (d xy and d x 2_ y 2)22i. In fact, we do observe this 
band weakly at other photon energies, like in Fig. (7^c), 
evidencing that it is not forbidden, but only weak and 
masked by the intense shallow band. 

On the other hand, the d xy * band, very clear in (al), 
is expected to be weak. As discussed in section II, its 
residual intensity is expected to be proportional to its 
hybridization with orbital antisymmetric with respect to 
z. In Fig. |6|a), we report the weight of d xz and d xy for 
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Figure 4. (a) Sketch of the Fermi Surface expected at k z —\. 
Thick and thin lines indicate main and folded bands in the 
sense of the lFe BZ with glide-mirror symmetry. Solid and 
dotted lines represent the in-phase and out-of-phase orbit als. 
Only in-phase orbitals are expected to have sizable intensity 
in ARPES (see section II). (b) FS measured at T=25 K, 
/kj=34eV and polarization along k y , except for the hole bands 
at Z, where it is in the xz plane, (c) 3D BZ for BaFe2As2, in- 
dicating the q wave vector for folding^, (d) Geometry of the 
experiment. The plane of light incidence is xz and the anal- 
yser slits are fixed along k y . (e) FS crossings around extracted 
from a (green) and ft (blue), (f) Unfolded Fermi Surfaces at 
k z — \ in lFe BZ with d yz spectral weight, (g) Same as (f) 
with all Fe 3d orbitals. 



this band, as obtained from our band calculation using 
Wien2k. We see that the d xz weight is in fact quite sig- 
nificant near k z =l, which may explain why it is clearly 
seen. Interestingly, it vanishes for lower k z value, giv- 
ing a way to check that the residual intensity indeed de- 
pends on this d xz weight. In Fig. [6^b), we follow the 
weight of this band at the Fermi level as a function of k z , 
by changing the photon energy. The band diameter de- 
creases away from k z , as will be analyzed in III-C, but it 
also sharply loses intensity and disappears in k z =0. This 
confirms very nicely the origin of this weight. The fact 
that it is a shadow band explains that it is not seen in 
most ARPES data. 

Finally, it is then a combination of spectral weight dis- 
tribution and matrix element effects that allows to isolate 
one electron ellipse. The very simple spectra allows to 
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Figure 5. (a-c) Dispersion observed and calculated along the 3 
cuts a, b, c, indicated in Fig. [4^e). The first row corresponds 
to panel /3 rotated by 90°, i.e. allowed bands are even/xz; the 
second row to a, i.e. allowed bands are odd/xz; the third row 
to DFT calculations; the fourth and fifth rows to unfolded 
spectral function in lFe BZ with only d xz and all d orbitals 
respectively. For (b4), because there is no d xz intensity in the 
original cut b, the spectral weight is plotted along a parallel 
k path 0.04 x ^ closer to Z and enhanced by a factor of 20. 
For the DFT calculation (third row), the style of the lines 
are the same as in Fig. 3(a). We show colors of the main 
orbital character for clarity, although there are non-negligible 
hybridization effects. White arrows in (cl) and (c2) indicate 
kinks at the band crossings. 



observe further the structure of the ellipse and how the 
bands couple together to form it. To understand how 
one goes from the deep electron band to the shallow one, 
it is instructive to consider the direction noted c, just 
next to the major axis. The hybridization is forbidden 
along k x because the two bands have opposite symme- 
tries, but becomes allowed along k y . Consequently, they 
cross without hybridizing in (a3), but, as soon as one gets 
away from this axis, as in (c3), a hybridization gap opens 
at their crossings. This forms an upper band, which pro- 
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Figure 6. (a) Weight for d xy and d xz in the deep electron band 
as a function of k z obtained in the calculation, (b) Intensity of 
the even deep electron band along ZX integrated at the Fermi 
level, as a function of k z . (c) Unfolded spectral function for 
the deep electron band shown in (b). 

gressively evolves into the shallow band of (b3). This 
is closely supported by our data where "kinks" signaling 
these hybridizations are clearly visible in (cl) and (c2) 
(white arrows), at the position of the two bands cross- 
ing, although we only image one band at a time. 

B. One-Fe perspective: unfolded band structure 

All our analysis so far are based on the 2Fe perspective, 
in which the phase interference between two Fe atoms in 
the unit cell are shown crucial to the observed ARPES 
spectral weight. Note however that in addition to the 
2Fe basis, the same information of quasi-particle exci- 
tation can be represented in two other basis sets: the 
lFe basis with glide-mirror symmetry, and the lFe basis 
with plain translation. These three basis sets correspond 
to Fig. [TJc), (d) and (e) respectively, which make obvi- 
ous the equivalence of the amount of information they 
host: compared to the 2Fe picture, the lFe pictures have 
only "half" of the bands, but reside in a twice larger 
k-space. The use of glide-mirror symmetry gives clean 
bands, while the use of plain translation gives more "par- 
tial" bands. Figure [T] should also make it obvious that all 
three representations respect the same sum rule and have 
the same total spectral weight. It is thus interesting to 
also consider the ARPES spectrum in an lFe perspective. 

Specifically, the lFe picture with plain translation has 
similar translational symmetry to the outgoing photo- 
electron, and would thus give a direct comparison to the 
ARPES spectrum. We thus proceed by unfolding the 
electronic structure in the lFe basis with plain trans- 
lation. Since the basis consists of a single Fe atom in 
the unit cell, the unfolded band structure absorbs the 
structure related phase interference, leaving only a triv- 
ial orbital dependent dipole matrix element of a single Fe 
atom, tunable by the polarization of the incident light. 

Figure |4jg) shows the calculated unfolded Fermi sur- 
faces in the lFe BZ, colored according to the orbital 
character and with intensity proportional to the calcu- 
lated spectral weight. This result is obviously very dif- 
ferent from the lFe picture under the glide- mirror sym- 
metry (c.f. Fig.JTJd)), but is in good agreement with the 



expected spectral distribution in Fig. [TJe). Given that 
our experimental polarization is along the ^/-direction, 
the atomic matrix element would thus be dominated by 
the yz orbital. Indeed, Fig. |4^f) gives the Fermi surface 
weighted by only the yz character and appears to repro- 
duce very well the experimental observation in Fig.[4^b). 

The power and simplicity of the unfolding picture is 
further illustrated in the last two rows of Fig. [5j which 
follow the similar /c-paths as those of the first two rows. 
Let's first focus on the last row (a5-c5) and only pay at- 
tention to the yz (green) character that contains very 
strong spectral weight. Evidently, all our experimentally 
observed features in the 2nd row (a2-c2) are in excellent 
agreement with the unfolded band structure. Interest- 
ingly, the calculation also establishes that the ARPES 
feature observed in (c2) is actually not from a single 
band, but from two crossing bands, as hinted by the ob- 
served "break" in the dispersion [indicated by the arrows 
in (c2)]. A similarly excellent agreement is found be- 
tween the fourth row (a4-c4) and the first row (al-cl). 
Particularly, the same "break" in the dispersion [indicated 
by arrows in (cl)] can be resolved from the two crossing 
bands. [For clarity and to be consistent with the analysis 
in the previous section, (a4-c4) shows only the weak in- 
tensity of the xz character on the electron pocket around 
(7r,0,7r), which gives the same intensity as the experimen- 
tally observed yz character around (7r,0,7r).] 

Notice that exactly along the path corresponding to 
Fig. [5|bl), there is practically no theoretical spectral 
weight from any of the even unfolded bands. In our ex- 
perimental spectrum, however, a very weak band can be 
recognized. A careful analysis of the unfolded band struc- 
ture indicate that this very weak feature most likely orig- 
inates from a finite momentum resolution of ARPES. In- 
deed, Fig. |5|b4) shows that slightly away from the path, 
a weak spectral weight appears that resemble the exper- 
imental observation. (An energy-broadening of the theo- 
retical spectrum is necessary to make a proper compari- 
son.) 

It is also instructive to compare the theory in the 2Fe 
picture in Fig. [5|a3-c3) and the lFe picture in (a5-c5). 
They confirm the previous analysis, as all in-phase bands 
have strong weight and the out-of-phase bands no or 
weaker weight with a different color. Of particular in- 
terest is the character of the electron pockets in (a3) 
and (a5). As shown in a recent study 18 , the reduced 
spectral weight at the Fermi level in this direction re- 
flects the important fact that the electron pockets are 
actually created by the breaking of the lFe translational 
symmetry. The weak (red and blue) electron pockets 
in Fig. [5|a5) appear to have the opposite character as 
the regular assignment in the 2Fe picture in (a3). For 
example, the lighter xy (blue) band in the 2Fe picture 
appears to be heavier than the xz (red) band in the lFe 
picture. Since the ARPES spectral weight will coincide 
with the unfolded band structure in (a5), an experimen- 
tal analysis of the polarization dependence might mistak- 
enly assign the heavier pocket as the xy pocket in the 2Fe 
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picture, seemingly in contradiction with the theoretical 
results in (a3). Based on such comparison, one might 
even draw the wrong conclusion of a strongly orbital de- 
pendent mass renormalization. However, this switch of 
character is exactly to be expected from the unfolding 
theory. As recently shown 18 , the folding potential here 
originates from the alternating out-of-plane positioning 
of As atoms, and thus only couples orbitals with op- 
posite z-parity. Thus, the folding of band structures is 
always "parity switching". Therefore, a xz pocket near 
(0,7r,0) would be folded to (7r,0,7r) and switch to xy char- 
acter instead. (The same "parity switch" would apply to 
the other xy pocket as well to give the folded pocket xz 
character.) This parity switching characteristic has to 
be taken into account when assigning band character in 
ARPES analysis, and highlights the convenience of the 
theoretical unfolded band structure in assisting ARPES 
analysis. 

Finally, not surprisingly, the k z dependence of the spec- 
tral weight can also be captured. Figure |6jc) shows that 
the growing xz character near k z = 7r, and consequently 
a strongly enhanced ARPES spectral weight, in perfect 
agreement with the experiment in Figure |6|b). It should 
be obvious that the unfolded band structure and Fermi 
surfaces offer a simple and direct way to analyze ARPES 
spectral function in general. It is particularly valuable for 
the case of families of Fe-based superconductors, due to 
the strongly k-dependent orbital character of the bands, 
and the special "parity switching" nature of the band fold- 



C. k z dispersion 

To completely characterize the band structure, we 
should follow the dispersion as a function of k z . The- 
oretically, the ellipse major axis is expected to rotate by 
90° between k z = 1 and k z = 0. This is required to have 
a 42 /m axis at X. This is obtained by squeezing the pock- 
ets, as sketched in Fig. [7^a) for one pocketful Note that 
the orbital character does not rotate, so that the asym- 
metry in vp between d xy and d xz /d yz remains the same 
at both k z . As the folded pocket has opposite characters 
(d xz along k x and d xy along k y ), the superposition of the 
2 pockets is needed to obtain the equivalency between 
the 2 perpendicular directions. 

Fig. [6^b) shows that the deep electron band diameter 
indeed decreases away from k z = 1 and Fig. [7^c) and 
(d) show that the shallow band diameter significantly 
expands between photon energies corresponding to k z =l 
and 0. In Fig. (7^b), we report by points fcf as a func- 
tion of k z for the deep electron band (measured along k x ) 
and the shallow one (measured along k y ). We compare 
these points to those expected in the calculation (thin 
lines). Qualitatively, the expected evolution is there, but 
it is quantitatively much weaker. In fact, both bands 
form significantly smaller pockets at all k z than the cal- 
culated ones; kp is nearly a factor 2 smaller. This ten- 




k -M1'A) k-k x (1/A) 

Figure 7. (a) Evolution of the electron pocket evolution ex- 
pected as a function of kz. (b) Size of the electron pockets 
as a function of k z . Points indicate /cf along k x for the deep 
electron band (blue points) and along k y for the shallow one 
(green points). The points for the deep band were extracted 
from Fig. 5(c). For the shallow electron band, two samples 
were measured for 20</kj<45eV (circles) and 40<hu<90eV 
(stars). Thick lines correspond to values expected in calcu- 
lations for BaFe2As2 doped by 10% Co. Thin lines are the 
contours calculated at -0.15eV. (c) Shallow electron band at 
39eV (k z [X]=7) (d) Shallow electron band at 27eV (k z [X]=6). 



dency was already noted both by de Haas- Van Alphen 
measurements in other compounds^ and ARPESPST We 
observe that the contours calculated at -0.15eV for both 
bands (thick lines) are in much better agreement with 
our data, supporting the idea that such a very large shift 
is necessary to adjust the data. It is very important to 
define this shift accurately to estimate the renormaliza- 
tion of the bands. With this shift, the warping expected 
along k z corresponds quite well to the calculated one. 



IV. CONCLUSION 

Finally, we obtained very simple spectra for the elec- 
tron pockets, despite the complexity of the band struc- 
ture. Only one band is detected in each measurement, 
d yz in the odd configuration and d xy * in the even one. 
These two bands combine to form one full ellipse in the 
sense of the lFe BZ of Fig. 1(d), as sketched in Fig. 3(e). 
To our knowledge, this is the first time that a complete el- 
lipse can be clearly resolved by ARPES. We explain why 
the other ellipse, expected to be formed by d xz * /d xy , is 
largely suppressed in our measurements by a combina- 
tion of spectral weight distribution (suppressing strongly 
d xz *) and matrix element effects (suppressing very effi- 
ciently d xy ). To clarify this point, we consider explicitly 
the spectral weight distribution expected from interfer- 
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ences between the 2 Fe of the unit cell and the associated 
parity switching effects. These effects were neglected in 
previous ARPES studies, but they are found very impor- 
tant to correctly interpret spectra near X points. Our 
data are in excellent agreement both with our qualita- 
tive symmetry arguments and with direct comparison 
with the unfolded band structure of refP^H We note that 
some authors prefer to describe electron pockets as one 
outer sheet (mainly of d xy character) and one inner sheet 
(mainly of d xz j d yz character) rather than as 2 crossed 
ellipses. The two descriptions can be consistent with our 
analysis, provided the gap between the two sheets re- 
mains quite small, as our measurements rule out a strong 
deviation of the electron pockets from the elliptic shape. 
We emphasize that, in any case, the spectral weight will 
be strongly modulated around these sheets, which will 
affect all electronic properties and disrupts the notion of 
sheet. 

Resolving one full ellipse is a favorable experimental 
situation, as we can then study in details the proper- 
ties of each bands, how they couple together and how 
they evolve with k z . As the other ellipse can be ob- 
tained by folding, there is no information missing about 
the electronic structure of the electron pockets. We con- 
firm unambiguously the orbital characters predicted by 
the band structure calculations for the electron pockets 
and a large anisotropy in Fermi velocities, the d xy band 
having a much steeper dispersion than d xz j d yz . We 
detect a k z dependence of the pockets qualitatively con- 
sistent with theory, although their sizes itself is much 
smaller. The calculated bands must be shifted by 0.15eV 
to match the experiment. 

More generally, our study allows to get a better under- 
standing of how to treat folding effects in ARPES. This 
should be useful for further studies of pnictides and in 
all, rather common, situations where the unit cell con- 



tains 2 inequivalent atoms or more. We show that the 
spectral weight distribution does not follow the 2Fe BZ 
schem e, nor that commonly referred to as the lFe BZ 
scheme^ 13 * 14 *. The key to understand the folding here is 
to consider separately bands having the symmetry of the 
folding potential or the opposite symmetry. An impor- 
tant consequence we derive and confirm experimentally 
is that the symmetry of the bands are not on-site prop- 
erties anymore, but that they depend on the location in 
k-space and on the nature of the band (main or folded). 

Furthermore, the distribution of spectral weight is an 
intrinsic feature of the electronic structure that defines 
its symmetry and can have important consequences. For 
example, the FS possesses a 4- fold rotation axis at the X 
point in the 2Fe BZ description (Fig. lc, this becomes a 
42 /m screw axis in the 3D structure of BaFe2As2), but 
only a 2-fold rotation axis in the lFe BZ (Fig. Id). Which 
is the correct symmetry in iron pnictides ? Recently, the 
anisotropy of spin fluctuations measured by neutrons at 
the X point were attributed to the true symmetry of the 
underlying Fe lattice that misses the 4-fold axiiPP. To 
what extent this "hidden" anisotropy can play a role in 
the various "anisotropics" and "nemacities" observed in 
these materials is still an open question. Our ARPES 
study clearly evidence the spectral weight modulation 
corresponding to the folding, implying there is no true 
42 /m axis at the X point. 
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